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Abstract 

A toy model is proposed which incorporates the reversible mode coupling mechanism re- 
sponsible for ergodic-nonergodic transition with trivial Hamiltonian in the mode coupling 
theory (MCT) of structural glass transition. The model can be analyzed without relying on 
uncontrolled approximations inevitable in the current MCT. The strength of hopping pro- 
cesses can be easily tuned and the ideal glass transition is reproduced only in a certain range 
of the strength. On the basis of the analyses of our model we discuss about a sharp ergodic- 
nonergodic transition and its smearing out by "hopping". 
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The mode coupling theory(MCT) as applied to supercooled liquid and glass transition is the 
only existing first principle theory for this last stronghold of condensed matter physics with sur- 
prising success in describing the initial stages of freezing|l|. However, as one expects, there are lots 
of formidable issues in applying the theory originally designed for critical phenomena dealing with 
thousands of Angstroms to the short scale problems of at most tens of Angstroms. In the origi- 
nal derivation of the idealized MCT self-consistent equation for the density-density correlator, the 
factorization approximation of replacing the four-body time correlation function by the product 
of two-body time correlation functions was introduced, which is totally uncontrolled and whose 
region of validity is unknown. Also nature of the non-ergodic state after the freezing transition is 
far from clear. The difficulty is compounded if one goes over to the so-called extended MCT where 
the rapidly varying momentum variable is introduced in the MCT scheme, which was necessary 
to partially cover the thermally activated processes important after freezing transition has taken 
place. However, such rapidly varying momentum variable can never be sensibly treated by MCT|^. 
Furthermore, physical picture of the hopping processes is still lacking. 

In the recent years, the ideas and the methodology developed in the spin glass community are 
brought in to deal with structural glass problems, which was pioneered by Kirkpatrick, Wolynes and 
ThirumalaiH]. Here mean field type toy models like p-spin models are being analyzed producing 
deep insights. However, in all these spin glass models the glass transitions are driven by non-trivial 
Hamiltonians (or free energy functionals) H . 

On the other hand, in the above-mentioned MCT for structural glass, the transition is driven by 
the reversible mode coupling mechanism, and does not require a nontrivial Hamiltonian. Consid- 
ering these circumstances we feel it important to develop a mean field type toy model that mimics 
this reversible mode coupling mechanism with trivial Hamiltonian and that can be easily analyzed. 
Proposition of such a toy model is the purpose of this communication. The model can be exactly 
solved and yield the self-consistent equation for the relevant correlator of the type familiar in the 



mode coupling theories of supercooled liquid and glass transition, where the strength of hopping 
processes can be readily tuned. 

Our model contains the TV-component density variable a^ with i — 1,2,3,- ■■ N and the M- 
component velocity field variable ba with a — 1,2,3,- ■ ■ M. In the end N and M will tend to 
oo keeping the parameter S* = M/N finite. The model equations of motion are expressed by 
the following Langevin-type equations for nonlinear random coupled oscillators with damping and 
noise: 

UJ 



ba = --yba-uj'^Kjaaj -=Jija{^^aiaj-T6ij) + fa (lb) 

V jV 

<fa{t)fp{t')> = 2-lT5al35{t-t') (IC) 

where the upper dot is the time derivative and the angular bracket is the average over the Gaussian 
white thermal noise denoted by fa ■ Here and after we use the summation convention over repeated 
indices, lo, 7 and T are the positive constant parameters of the model and Ki^ are the coefficients 
satisfying the conditions KiaKip — Sap ( but not the conditions KiaKja — Sij for M < N). Jija 
are the mode coupling coefficients which are chosen to be "quenched" Gaussian random variables 
with the following properties: 
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where ^"^"^ denotes average over the J's. 

(Ilk) is analogous to the equation of continuity of fluid and (Elb) is like the equation of motion 
where the rhs is like the force acting on a fluid element. The coefficients Kia govern linearized 
reversible dynamics of the oscillators with the dynamical matrix ft given by flij = uj^KiaKja- 
In constructing this model we were motivated by |5|, |6| where random coupling models involving 
infinite component order parameter have been shown to be analyzed exactly by mean field type 
concepts. 

One can derive from the Langevin equations IW the corresponding Fokker-Planck equation for 
the probability distribution function Dt{{a},{b}) for our variable set denoted as {a}, {5}, whose 
stationary (equilibrium) solution is given by 

De{{a},{b}) = cst.e ^.=1=^"^ Z.c=i ^t". (3) 

where est. is the normalization factor. 

For the subsequent analysis it is most convenient to adapt the generating functional method 
of B to our case, which produces no new problem and its details can be omitted. Then below we 
obtain for this toy model the correlation functions defined by 

Cait,t') = —<aj{t)a,{t')>, Cab{t,t') = —Kja<aj{t)ba{t')>, (4a) 

CUi.t') = j^Kjo. < ba{t)a,{t') >, Cb{t,t') = ^ < ba{t)ba{t') > (4b) 



C^it,t') ^ -<a^W«a(0>, ci^^K.aa, (4c) 
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The last one is needed to close the self-consistent set of equations for correlators when M < Nm. 
We can write down the set of effective linear Langevin equations resulting from the generating 
functional method which are valid in the limit of M, TV — > 00: 

d^it) ^ K,aba{t)-^aai^a,{t)^K,a^abi^ba{t)+ft{t) (5a) 

bait) = -jba{t)^LU^K,aa^{t)~K,a^ba6^a,(t)-^bb6dba{t)+fi{t) (5b) 



where /" and /^ are the effective thermal noises whose correlations are related to the S's in a 
usual way, and 
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T,Q^a{t)= dt'T,(t,t')a(t') etc. (6) 

'^ — OO 

Here the kernels H's are given by 

Eaa{t,t') = ^S*{Ca{t,t')a{t,t')+5*Cab{t,t')CUt,t')), l^abit, t') = -2^6* Ca(t, t')CUt, t') 
Eba{t,t') = -2^S*Cait,t')Cab{t,t'), S,, (i, t') = ^^C, (i, t'f . (7) 

Using this effective Langevin equation we can derive the equation for correlation functions (y) 
which depend only on i — t', or their Laplace transforms the C^(z)'s (the S^(z)'s are the similarly 
defined Laplace transforms of S(t, t')'s): 

^Cf(z) = ^ + (l-^^,,(z))d*Ct{z)-^L{z)Ct{z) (8a) 

zCUz) = -il + ^idz))CUz)-{io^ + ^Liz))C^''{z) (8b) 

zCii(z) = (l-S^,(z))C,^(z)-Ei(z)Ci;(z) (8c) 

zC,^(z) = T-{cj' + j:ti^))C!:,{z)-{j + j:^,{z))Ct{z) (8d) 

zC^Hz) = ^ + {l-^^,b{z))Ct{z)-^L{z)C^\z) (8e) 

Since the "self-energies", the S's, are expressed in terms of the correlators by (|^), the equations 
^ constitute the self-consistent scheme for the 5 correlators (Q). These have to be solved numer- 
ically under the appropriate initial conditions for the C(f',t')'s. Eliminating C^(z) and C^^{z) 
from (m , b and e) we find the following equation for C^{z): 



Ca{z) = -2 



^2z + Ei„(z) 



1 - uj'5 



{^-Kbi-AY 



[z + ^L{z)){z -f 7 + S..(z)) + c.2(l - ^L^[z)y 



(9) 



where we have used the symmetry relation S^^(z) = —uj^Y.^^^{z), and other correlators entering 
the E's have been expressed similarly. 

For (5* = the {a} variables are time-independent since there is no velocity variable {6} that 
drives dynamics of {a}. Therefore the system is trivially nonergodic: Ca{t^ t') = Ca{t' , t') = TjuP'. 
For (5* = 1 (^) reduces to the equation derived in [|o| and the system is ergodic due to the presence 
of strong hopping process. For intermediate values of (5* we may expect ergodic-to-nonergodic 
transition as we have numerically seen for 8* — 0.3 (see Fig.l). For small values of 8* we have the 
following result correct up to first order in 8*: 
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1-5*'+^'^ 



(10) 



1 + 2g^T 
The phase diagram in the variables 6* and T found numerically is displayed in Fig. 2. Equating 



(10) to zero we find the transition at T^ = 5^'^(1 — S*)/{6* — 2) < 0, the unphysical value. Hence 
we suppose that Tc should diverge at some nonvanishing positive value of 6* . 

It should be emphasized that our theory is exact in the limit of large M and N. By varying 
the parameter S* which is the ratio of the numbers of components of {a} and {6}, we can control 
the strength of the hopping processes represented by the T,aa and T,ab terms. Reducing S* drives 
the system into glassy region. 

Our naive expectation when we started this work which was for S* = 1 {i.e. M = N and 
Kid = Sia) was the following. We are constructing an TV-component model for which mean field 
mode-coupling approximation is exact as N —^ oo. A mean field type approximation often gives 
a sharp transition and Gotze' s MCT without "hopping" can be regarded as a kind of mean field 
theory. Hence it is natural to expect that this toy model designed to reproduce the Gotze type 



mean field theory rigorously has such a sharp transition. But unexpectedly our toy model did not 
show a sharp transition. In order to produce a transition it was necessary to generalize the original 
A^-component model to allow M < N . Thus, understanding this so-called ideal glass transition 
without relying on uncontrolled approximation appears to be more subtle than expected. Another 
related issue is concerned with the nature of the "hopping" process which can occur even with 
trivial Hamiltonian. This means that the "hopping" found in our model has nothing to do with 
thermally activated energy barrier crossing since the trivial Hamiltonian does not possess such a 
barrier. The same comment applies to [|lO| . Thus we are faced with the problem of understanding 
the Arrhenius type hopping rate found recently |ll|. 

One may interpret that since the reduction of 5* is tantamount to restricting accessible region 
of the phase space, we can say that we have an entropic barrier [ll2| . This interpretation, however, 
need to be substantiated. In this connection it is illuminating to eliminate the variables {6} 
adiabatically to obtain a closed Fokker-Planck type equation containing only the {a} variables. 
Then we can show that the diffusion marix in this equation is singular |13[, i.e. its determinant 
vanishes. This gives rise to a possibilty that the equation can have stationary solutions other than 
the equilibrium Gaussian oneQ. Such solutions are precisely the kind of non-ergodic states found 
numerically for the present original model. This feature will be signifincant beyond the present 
toy model. We also note the simultaneous appearance of the ergodic to non-ergodic transition and 
the new variable a^ = KjaUj, (^). Elucidating its implication is an intriguing future task. 

We hope to discuss these aspects further including critical behavior near transition and an 
extension to aging in future. 
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Figure Captions 
Fig.l Temporal behavior of Ca{t) = Ca{t,0) ioi g = j = uj = 1, S* = O.3.. The curves are, from 
left to right at long times, for T = 1.0, 0.5, 0.3, 0.2, 0.1, 0.05, 0.01, 0.001. 
Fig. 2 The phase diagram in the S* — T space dividing the ergodic and non-ergodic regions. 



g=^co=l, 5*=0.3, 1=1.0,0.5,0.3,0.2,0.1,0.05,0.01,0.001 



0.8 



0.6 



U 



0.4 



0.2 












10 



100 



1000 



10000 



0.06 



0.04 



ERGODIC PHASE 



0.02 



NONERGODIC PHASE 







0.2 0.4 0.6 0.8 

5* 



